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Abstract 

Recent development in exploring flavour dynamics in the supersymmetric extension of the Stan- 
dard Model is reviewed. Emphasis is put on possible interesting effects in 6-physics arising for large 
values of tan (3 both in the case of minimal flavour violation and in the case of flavour violation origi- 
nating in the sfermion sector. The importance of the flavour changing neutral Higgs boson couplings 
generated by the scalar penguin diagrams and their role in the interplay of neutral _B-meson mixing 
and ^ M^M" decays is discussed. It is pointed out that observation of the fJ-'^ decay 

with BR at the level 3 x 10~* would be a strong indication of nonminimal flavour violation in the 
quark sector. Possible impact of flavour violation in the slepton sector on neutrino physics is also 
discussed. 



^Talk given by P.H.Ch. at the International Meeting Physics from the Planck Scale to the Electroweak Scale, May 2002, 
Kazimierz Dolny, Poland and dedicated to Stefan Pokorski on his 60*'' birthday. To apper in the special volume of Acta 
Physica Polonica B. 



1 Introduction 

Physics of flavour and of CP violation continues to be an interesting subject to study in various extensions 
of the Standard Model (SM). On one hand, before the advent of LHC and linear colliders, which will 
enable us to probe energies much above the electroweak scale directly, rare processes intensively studied 
in numerous experiments are the first place where the effects of new physics - i.e. virtual effects of new 
particles - can most likely be detected. On the other hand, studies of baryogenesis [0 strongly suggest 
that the SM with its unique source of CP violation and known particle content is unable to explain 
the baryon to photon density number ratio, n^/n^ 10~^, observed in the Universe, given the present 
lower limit on the Higgs boson mass. This supports expectations that some deviations from the SM 
predictions will eventually be encountered in the ongoing or planned precision studies of rare and CP 
violating processes. Finally, the first direct indication of inadequacy of the SM has also to do with flavour 
physics - namely with neutrino oscillations. Explanation of the observed neutrino oscillations requires 
the introduction of flavour mixing (and perhaps also of CP violation) in the lepton sector. While the SM 
can easily be extended to describe neutrino oscillations, there are strong theoretical arguments that the 
observed phenomena have their origin in physics at energy scales much higher than the electroweak scale. 

Physics of flavour in the framework of supersymmetric extension of the Standard Model was also 
always in the center of Stefan Pokorski's interest. It is therefore a pleasure to devote this article to him. 
The subject is of course too vast to be reviewed here in all details. Instead we concentrate on its most 
interesting in our opinion aspects. These include recent investigations of the supersymmetric effects in 
&-physics arising for large tan /3 and in the neutrino sector. In this context it is appropriate to recall here 
that systematic investigations of supersymmetry at large tan/3 begun with the Stefan Pokorski's seminal 



2 Flavour violation: minimal and generalized minimal 

Extensions of the SM can be divided into two broad classes: models in which the Cabibbo-Kobayashi- 
Maskawa matrix (CKM) in the quark sector and Maki-Nakagawa-Sakata matrix (MNS) in the lepton 
sector are the only sources of flavour and CP violation and models in which there are entirely new 
sources of flavour and/or CP violation. Both options can be realized independently in the quark and 
lepton sectors of the simplest supersymmetric extension of the SM - the MSSM - and it is the experimental 
task of utmost importance to establish which one is realized in Nature. 

If the CKM matrix is the only source of flavour and CP violation in the quark sector, the natural 
question is what is the impact of new physics on the determination of its elements. In particular one wants 
to know the value of the Vtd element which is needed e.g. to predict the rate of the decay — > fi^fj,~ 
and other interesting rare processes. It is also important to see if the consistency of the determination of 
the CKM matrix elements from different processes imposes any constraints on the MSSM parameters. 

The CKM matrix V is most conveniently parameterized as follows M : 



The Wolfenstein parameters A « 0.222 ± 0.0018 and A « 0.83 ± 0.06 are rather accurately determined 
from transitions dominated by tree level contributions, and are hence insensitive to new physics. At 
present, processes of this kind put also some constraints on the remaining two (conveniently rescaled 
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Q) parameters p = p{l — and 77 = 77(1 — A^/2). The value of the combination i?b = \/ p"^ + rf 

is constrained to 0.27 < Rb ^ 0.46 by the result |K6|/|Vcfc| = 0.08 ± 0.02 (at 95 % C.L.) extracted 
from the charmless B decays. The CP violating time dependent asymmetry in the B — > ^Ks decay 
constrains the phase /3ut of the Vtd element: Vtd = |Vtd|e'*''"' . In minimal models this asymmetry is 
simply given by sin2/3ut. The average of the measurements done at BaBar and Belle gives sin2/3ut = 
277(1 — p)/ (1 — PY + = 0.78 ± 0.08 [||. There are also prospects for extracting from such processes 
also the phase 7ut of the Vub element: Vub ~ iKtble However, large theoretical and experimental 

uncertainties still prevent precise determination of p and 77 exclusively from tree level dominated processes. 

Parameters p and 77 are also extracted from measurements of the i?^ ^.'^1 d meson mass differences 
AA/g^d and of the parameter ek of CP violation in the neutral kaon system. This allows to overconstrain 
the values of p and fj and test the assumption of minimal flavour and CP violation in the quark sector. 
However, since all the three quantities are loop induced, the new physics can contribute to relevant 
amplitudes. Therefore, the values of p and 77 determined from AM^-^d and £k can significantly depend 
on new physics. One can also expect that consistency of the CKM parameters determination puts some 
constraints on new physics. 

To compute AMg^d and ek one integrates out from the theory all the states with masses }t Mw and 
constructs the effective Hamiltonian of the form 

T^eS = [q^^ ^CKM^xOx (2.1) 

where Ox are the local four-quark operators {X labels different Lorentz structures: X =VLL, VRR, 
VLR, SLL, SRR, SLR, TL, and TR) and ^ckm^^ ^'"^ their Wilson coefficients. An important feature 
of the minimal flavour violation is the factorization of the Wilson coefficients into A^p^f^ which depends 
only on the CKM matrix elements and Cx which to (a good approximation) are real numbers. 

At the level of the effective Hamiltonian (|]^) models of new physics in which the CKM matrix is the 
only source of flavour and CP violation in the quark sector can be further divided into two broad classes 
§: 

• the MFV (minimal flavour violation) models - truly minimal ones, in which, just as in the SM, only 
CvLL Wilson coefficient is non-negligible and Cvll responsible for ^^'^ K^-K" meson 
mixing are all equal (universal value of Cvll)- 

• the GMFV (generalized minimal flavour violation) models - in which more Cx are non-negligible 
and/or are non-universal. 

As we shall see, the MSSM can be of either type, depending on the ratio V2/V1 = tan/3 of the vacuum 
expectation values of the two Higgs boson doublets. 

Basic formulae used to determine p and 77 read (see e.g. refs. ^ ^ for further details): 

77 [(1 - p)A^V2F' + Pc] A^Bk = 0.204 , (2.2) 

where the number on the rhs stems from the measured value Ek = 2.28 x 10"'^, and 

AA/, = -^-^Mb,vbBb,FI^\V:,Vu?\F'\ « Bs,F|J(1 - p) - zn\^\F'\ 

f<2 jur2 

AMs = ^-^Mb^iibBbFI \V;,Vts?\F'\ oc BbFI . (2.3) 
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In eqs. (2^), (2_^) ?/2 — 0.57 and j]b — 0.55 summarize the short distance QCD corrections to Cvll 
Wilson coefBcients and Pc = 0.30 ± 0.05 is the known charmed quark loop contribution to ek- Factors 
Bk w 0.85±0.15, Bb^FI^ w (230±40 MeV)^ and Bb^F^^ w (265±40 MeV)^ § parameterizing matrix 
elements of the standard VLL operators are the biggest sources of uncertainties. The three factors F^, 
F'^ and F" can be expressed in terms of the Wilson coefficients Cx, their QCD RG running and matrix 
elements of the operators Ox for X ^ VLL. In a concrete model of new physics such as e.g. the MSSM, 
F^ , F"^ and F'* are calculable functions of its parameters. [| The distinction between MFV and GMFV 
models is reflected in that in the formers = F** F'^ whereas in the latter models all can be 
different. In the SM ^ F'^ ^ F" ^ Fsm = 50(7714) « 2.38 ± 0.11 for mt{mt) = 166 ± 5 GeV. 

The measured B^-B^^ mass difference, AMd = 0.496/ps, puts the constraint^ 1.04 < y^\W\Rt = 
•\/|F^|l — p—ifj\ ^ 1.69. The role of AM^ which does not depend directly on p and fj is twofold. Firstly, 
as follows from eqs. (2^), any new physics model must be such that F* it gives rise to satisfies |ll| 

F- 



< 



/ AMs \ 



(2.4) 



Since at present only the lower limit on the B'^-B'^ mass difference is known, AMg > 15/ps, the factor 
F'* is bounded only from below. Secondly, once measured, AM^ combined with AMd will allow for more 
precise determination of \Vtd\ oc \l ~ p — ifj\ because the ratio of F^ Bb^ to F^^Bb^ is known with 
better accuracy than these factors individually: ^ ~ 1.15 ± 0.06 ||l^. For given F'^ /F'\ the value of Rt is 
then determined, from the formula 



Rf. 



ps 



Note that Rt determined in this way is universal in the whole class of MFV models for which F'^ /F^^ = 1. 
In contrast, in GMFV models the extracted value of Rt does depend on new physics contributions to F** 
and/or F"^. 
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Figure 1: Contribution of the chargino-stop and charged Higgs (VF^) boson box diagrams to F^, F'' and 
F^ in the MSSM. Crossed diagrams are not shown. 



^AU the matrix elements of the operators Ox, also for X ^VLL, are now known from lattice calculations [y, |lO[ . 
Nevertheless, the uncertainties in their values still introduce some uncertainty in the factors which depend, apart from 
Wilson coefficients and calculable QCD RG factors, also on the ratios of these matrix elements to the matrix element of 
the standard VLL operator. 

■^AU bounds and allowed ranges of various quantities quoted in this article are obtained by scanning over all uncertainties 
within their respective Icr ranges. 
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3 Supersymmetric contributions to F^, and F' 



Dominant supersymmetric contributions to F'^ , F'^ and F"^ for small and moderate values of the ratio of 
the two MSSM Higgs boson doublets vacuum expectation values V2 /vi = tan /? are well studied |l^ . 
They arise from box diagrams shown in fig. |l| and for 2 < tan/3 < 20 give F^ = F'^ = F^ = F. Thus, for 
not too large values of tan/? the MSSM is of the MFV type. Maximal values F/Fsm ~ 1-4 are reached 
for lightest sparticles still not excluded by direct supersymmetry searches and tan /3 as small as possible.^ 
With increasing tan/3 and/or increasing sparticle and charged Higgs boson masses the value of F/Fsm 
decreases to 1 [|4[ p^ . 



(dfl).j (di)/ {dL)j {dR}! {dL)j {dR)i 




{dL)i {dR)j {dR)i {dL)j {dL)i [dR)j 



Figure 2: Double penguin diagram contributing to Cf^^^ Cf^^ and Wilson coefficients, respectively 
in the MSSM with large tan /?. 

As has been found recently |ll| (see also jl^), for large values of tan/3, ~ 50, the Wilson 
coefficient Csll , Csrr and Cslr of the effective Hamiltonian (^) can receive very large contributions 
from the so-called double scalar penguin diagrams (formally two-loop) shown in fig. ||. The origin of 
the fiavour changing couplings of the neutral Higgs bosons can be easiest understood in the effective 
Lagrangian approach |l^, |l^ (see also [Q): due to the triangle (scalar penguin) diagram shown in 
fig. ||a, integrating out sparticles (but not the Higgs bosons) in the approximation of unbroken electroweak 
symmetry generates the Yukawa coupling of the Higgs doublet to down-type quarks that is not present 
in the original MSSM Lagrangian. Thus, in the low energy effective Lagrangian both Higgs doublets, H'^ 
and i?", couple to down- type quarks and this, after the electroweak symmetry breaking, gives rise to the 
tree level flavour changing couplings of A", /i° and 

iff" 

A 

/ V q 

*■ ' ■* X *■ ■* 

dj g g II 
c) 

Figure 3: Diagrams generating flavour changing neutral Higgs boson couplings, g, u'^ and d'^ denote 
electroweak eigenstates. Diagrams b) and c) contribute only in the case of non-minimal flavour violation 
arising from squark mass matrices (sec. 5). 

^Recall also, that for the top squarks hghter than ^ 1 TeV the range 1 tan/3 2 is excluded by the unsuccessful 
search of the lightest Higgs boson at LEP. 
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For the transitions djdj ^ djdj the double penguin diagrams give 

R Qfcm JTif ,^2 A r, ( cos^ a sin^ a sin^ /3 

^•(^SRR jg obtained from C^ll replacing m^^ by m^^). is given by 
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where ' — M? /m^,, i,j — 1,2 are the ratios of the stop and chargino masses squared, the matrices 
Z+ and are defined in ref. and 

x\nx y^T^y 



H2{x,y) 



{l-x){x-y) {l-y){y-x) 



Because for Ma > Mz and tan/3 > 1 one has Mjj « M\, sin a w 0, the coefficients C^^^ and C^^^ 
are suppressed It turns out however, that for sufficiently large stop mixing parameter At the double 
penguin contribution to (7^'"^ for the hs sb transition is significant despite the suppression by the 
strange quark mass. Inserting numbers one finds 



tc 



for nib = 2.7 GeV, nis — 60 MeV at the scale Q — rrit, Mh — Ma and sin a — 0. For tan/3 ~ 50, Xtc ^ 
0(1) and CP-odd Higgs boson not too heavy this is comparable with the value of the Wilson coefficient 
of the standard VLL operator: C^ll ^ tiSoirfit) « 9.5. The ratio (^SLr^^^vll further increased by 
the QCD RG effects ||: C^^^{A£ GeV) = 2.23 C^^^{mt) while C^^^{A.Q GeV) = 0.84 C^^^(mt). For 
the transitions hd db and ds sd similar double penguin contributions to C'S'"^ are negligible being 
suppressed by uid/ms « 0.06 and md/rrn, « 0.001, respectively. Thus, for large values of tan/3 the MSSM 
becomes of the GMFV type with ^ F'^ Fsm + and F'/Fsm < 1- 

The important features of the double penguin contribution to C^^^ are the following: it grows as 
tan^ (3, it is always negative leading to F^ < Fsm and is directly sensitive to the top squarks mixing 
(-(^SLR (X ^j). Moreover it does not vanish if all the sparticle mass parameters are uniformly scaled up 
(non-decoupling effect). It does however vanish as the inverse square of the Higgs sector mass scale set 
by Ma. 

Figure ^ showing constraints from different experimental data in the (p, fj) plane allows to discuss the 
value of Vtd in the two scenarios: MSSM with small and large tan/3 as a function of measured in the 
future value of AAfg. {Rt oc \Vtd\ equals the length of the line connecting a given point in the (p, fj) plane 
with the point (1,0).) 

In the MFV-type MSSM with small tan/3, and also in the SM, F'^/F'^ = 1 and p and fj are bound 
to lie inside the black spots in figure || which are compatible (for sin2/3ut ~ 0.78) with the constraints 
imposed on Rt by the value of iKib/Kbl- Therefore, Vtd determined from |y„b/Vc6| (oc Rb), AAIg/ AAId 
(oc Rt) and the asymmetry measured in the B ipKs decay (= sin2/3ut) in the MFV-type MSSM and 
in the SM is the same: \Vtd\ = (7.75 - 9.5) x lO^^ {R^ = 0.90 - 0.99) for AM, = (15 ± 0.5)/ps and 
\Vtd\ = (6.7-8.2) X 10"^ {Rt = 0.78-0.85) for AM, = (20±0.5)/ps. Taking into account the constraint 
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Figure 4: Ranges of (p, rj) allowed at Icr for AMg = (15.0 ± 0.5)/ps (upper panel) and (20.0 ± 0.5)/ps 
(lower panel) for sin2/?ut = 0.78 ±0.08 and different values oi /F'^ (marked in the figures). Black spots 



correspond to F'^/F — 1. Dotted (dash dotted) lines show the constraint on (p, ?/) from ex (eq. (2.2)) 
for F'^ — FsM (for F'^ ~ 1.3Fsm in the upper and for = l.Si^sM in the lower panels, respectively) 
Solid semicircles mark the range of Ri, = \/p^rf)^ allowed by iKib/Vchl. 
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imposed on p and fj by Ek does not change anything for AM^ close to 20/ps (the shaded region hes 
entirely between the two Ek hyperbolae even for F"^ — 1.5). On the other hand, if the value of AA/^ is 
close to its present lower limit of 15/ps, it follows from ( |2.4| ) that — F'^ — F^ must be smaller than 
« 1.3 (MSSM parameters leading to i^^ bigger than 1.3 are excluded). ^ Only for close to 1.3 can the 
upper edge of allowed Rt values (and therefore of \Vtd\) determined from the fit to the data be slightly 
lower than in the SM. We conclude therefore that from the practical point of view the value of Vtd in 



the SM and in the MFV type MSSM is the same. Note also, that in this scenario the factors Fb^ \ B 



and Fb^ y a-re positively correlated in the sense that, for fixed AAfg, bigger values of F'^ — F'^ — 
require both these factors to assume simultaneously values from the lower parts of their respective ranges 
obtained from lattice calculations. 

In the MSSM with tan/5 ^ 50 F" ^ F'^ = Fsm and F'^/F'^ = F^/Fsm < 1- The absolute bound 
does not allow for \F'' /F'^\ < 0.5 for AAIg — 15/ps but the inspection of the upper panel of figure |4| shows 
that the combination of constraints imposed on p and fj by Ek (dotted lines) and Rf, (solid semicircles) 
excludes also those MSSM parameters for which I^VF^^I = |FVFsm| < 0.55. Similarly, for AM^ = 20/ps 
the bound dJ) gives \F'/F'^\ > 0.69 whereas Ek and Rt require \F'/F'^\ > 0.75. For values of \F''/F'^\ at 
the lower edge of the allowed range the value of |Vtd| extracted from AMg/ AMd is smaller than in the SM 
(for example, for AM^ = (15±0.5)/ps and IF'/F'^] = 0.6 or for AM^ = (20±0.5)/ps and \F'/F'^\ = 0.8 
\Vtd\ = (6.0 - 7.3) X 10"-''). Note however, that for -1 > F^/F"^ > -1.3 for AM^ = (15 ± 0.5)/ps 
(-1 > F^/F"^ > -1.76 for AM, = (20 ± 0.5)/ps) the value of \Vtd\ can be bigger than in the SM. Of 
course, large departures of |F''/F'^| from 1 discussed here are compatible with AAIg and AMd separately 



provided the lattice factors Fb^ \ Bb^ and Fb^ \ Bb^ assume appropriate values (which however remain 



within their respective uncertainties). The correlation of Fb^\/ Bb^ and Fb^\/ Bb^ is again positive 



(although weaker than in the previous case): smaller F'^/F'^ — F'^/Fsm requires bigger Fb^\/Bb, (to 



reproduce AM^) and leads to smaller value of 114^1 which in turn calls for bigger Fs^y Bb^ to reproduce 
AMd. 



4 Impact of the scalar penguins 

As has been pointed out in ref. |23j, the reliable calculation of the flavour changing neutral Higgs boson 
couplings in the MSSM requires resummation of the tan (3 enhanced terms. Furthermore, as demonstrated 
in refs. ||2j, |25| there are also tan/3 enhanced corrections to the couplings of the charged Higgs and 
Goldstone bosons which affect the box diagram contribution of these particles to the Wilson coefflcients 



of the effective Hamiltonian (2.1). Technical details and systematic study of all these refinements can be 
found in |^ and will be not discussed here. They are however included in the numerical results presented 
below. 

The role of the scalar penguin induced flavour changing neutral Higgs boson couplings is twofold. 
Firstly, for tan (3 > 30 a big portion of the MSSM parameter space (the bigger the higher is the lower 
experimental limit on AM^) in which the parameter At is large (and hence the mixing of left and right 



top squarks is substantial) is excluded by the bound (2.4) and its reflnement related to constraints on p 



and r] from ek and |T4b/Kb| discussed in the preceding section. Typical dependence of F^/Fsm on the 

*Note that this puts severe constraints on the scenario with tan/3 < 1: stops, charginos and /f+ would have to be very 
heavy in order their contribution to B^-B^ mixing described by be sufficiently suppressed. 



7 



MSSM parameters is shown in figure |5|. For /i > the resummation of tan j3 enhanced terms mentioned 
above increases the value of f /i^sM (i-e. suppresses the negative contribution of the flavour changing 
coupHngs of neutral Higgs bosons) compared to the naive one- loop calculation of ref. |2^. For /i < 0, 
however, the effects of the flavour changing couplings are enhanced by the resummation. The parameters 
in flgure ^ has been chosen so that ^At has always the sign which allows for cancellation of the 
and chargino-stop contributions to the amplitude of the B — > X^^ decay. 




10 20 30 40 50 ■ 200 400 600 800 1000 



iau§ Mh+ (GeV) 



Figure 5: F"^ /F^m as a function of tan (3 and Mu+ for the lighter chargino mass 750 GeV and \M2/ lj\ = 1- 
Solid and dashed dashed lines correspond to stop masses (in GeV) (500,850) whereas dotted and dot- 
dashed lines to (600, 750). The mixing angle between the two stops is \6f\ = 10°. Solid and dotted 
(dashed and dot-dashed) lines correspond to /i < (/i > 0) and the stop mixing angle has the sign 
opposite to that of /i. nig = 3M2 and the right sbottom mass is 800 GeV. 

Secondly, the same flavour changing neutral Higgs boson couplings which (through the double pen- 
guin diagrams) affect the B'^-B'^ mixing has been found |l9|, to totally dominate for tan /? > 30 
amplitudes of the decays B°^ ^^^fJ■^ ■ Calculating the diagram shown in fig. O one finds Hi 

A{B"g -> fi+fi-) = u(fci) (6 + a-f') v{k2) (4.5) 

where q — s or d, u(fci), v{k2) are spinors of the final state leptons, and (without resummations, with 
Mh « Ma etc.) 



Gf (Xcm ^'^B„ 



« = 5 = -V,,V,,.nF,^^^^X,ct.n p . 

Therefore, in the MSSM with large tan/? the decay rate behaves as BR{B'^^ fi^fi^) oc (tan^/3/M^) 
and - without additional constraint imposed - could, for tan /3 > 50 and the Higgs bosons not too heavy, 
even exceed the present experimental bounds |^ 

BR{B°^ fi+^r) < 2.0 X 10^^ CDF 

BR{B^ < 2.1 X 10~^ BaBar (4.6) 
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That is, the rates predicted in the MSSM could exceed by 3-4 orders of magnitude those of the SM 



BR{B'^ ^ /i+A*")sM ~ 3.5 X 10 
BR{B° ^ ^i+^I-)sM ~ 1-4 X 10-1" 



V 230 MeV 

Fb, 



200 MeV / V 0.009 




Figure 6: Flavour changing neutral Higgs boson couplings contribution to the amplitude of the B^ ^ ^ 
decays. 

However, as we have discussed above, for light A" the magnitude of the flavour changing scalar 
coupHngs bftA^SL and baH^SL (and, hence, also of the couplings b^A^dL and bnH^dL, because in the 



GMFV MSSM they are proportional to the former ones) is strongly constrained by the condition (2.4). 
Therefore one can expect that also the contribution of the neutral Higgs boson exchange shown in figure 
^ to the amplitudes of the B^ M^/^^ decays is bounded by the condition (^.4|). In other words, the 
lower limit on AM^ should put the upper bound on the possible values of BR{B1 ^ predicted 
in the MSSM. 

Figure ^ shows the correlation of the predicted values of BR{B^^^ ii^ji^) and AA/s for a sample 
of the MSSM parameters for Ma = 200 GcV and tan/3 50. In the case of the BR{B^ M^M") we 
have determined the value of \Vtd\ consistently, that is we have scanned over the Wolfenstein parameters 



A, A, p and fj as well as over the nonperturbative parameters Fb^ y Bb^ and computed the rate only for 
those A, A, p and 77 for which eif, AMd, sin2/3ut, IKifj/Kbl assumed acceptable values. We have also 
excluded all points for which the rate of the B Xgj is unacceptable. 

The upper bounds on BR{B^ ^ — > p,^p^) are clearly seen in figure |7[ For tan /3 = 50 and Ma = 200 
GeV all points for which F'/Fsm < -0.52 (so that AM^ > 15/ps) give BR{B'^^ M+M") above the CDF 
bound (4^) and excluding also points for which AMg < 15/ps we see, that BR{B^ — > p'^p,^) < 10^^ 



and BR{B^ p~^p ) < 3 x 10 ^. Points for which /Fsm < —0.52 can survive for smaller values of 
tan/3 and/or heavier CP-odd scalar A° (note that BR oc tan^ /?/M| whereas AM^ oc IF'*] cx tan'' /?/M|). 
In this case however both, AMs and BR{B'^^ p^p~) are entirely dominated by the contributions of 
the scalar penguins and it is easy to estimate that whenever BR{B^ p^p^) is below the CDF bound 



( M3|) , BR{B° p+p^) ;5 6 X 10"^ i.e. it is below the BaBar bound (gj). 

We conclude that the MSSM parameter space in which the parameter At is not unnaturally big (that 
is. At ^ Msusy) is more strongly constrained by the lower limit on AMs than by the non-observation 
of the B^^^ decays in CDF and BaBar. In particular the bound BR{B'^ — > M^A^ ) < 3 x 10"* 

holds. For parameters such that F'^ /Fsm < —0.52 (larger At) there is a weaker upper bound BR{B^ 
p+p-) > 6 X 10-* 
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Figure 7: Correlation of BR{B'^ fi^fj, ) (left panel) and BR{B^ fi^fj, ) (right panel) with AMs in 
the GMFV-type MSSM for tan/3 = 50 and Ma = 200 GeV. In this case all points for which FV^sm < 
—0.52 (so that AMg > 15/ps) give BR{B'^ fi^^~) above the CDF bound (4.6) and have been discarded. 
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5 Flavour violation in squark mass matrices 



In supersymmetric extension of the SM, flavour and CP violation can originate also in the sfermion sector. 
In general, the 6x6 mass squared matrices of left- and right-chiral sfermions of the same electric charge 
have the form^ 

m9j ) ( 



where {m'^^ etc. are 3x3 submatrices. If the latters are not diagonal in the so-called superCKM 
basis, in which quark mass matrices are diagonal, then their off-diagonal entries generate flavour changing 
neutral currents. For example, large, oc a^, contributions to K'^-K^ or B'^ ^-B'^ ^ mixing are then induced 

by the gluino box diagrams shown in figure ^. In this figure the off-diagonal entries of matrices (j^Ixy^ 
are treated as additional interactions (the so-called mass insertion approximation pOl ^ fH)- As these 



contributions are not proportional to the CKM matrix factors the effective Hamiltonian (2A) for |AF| — 2 
transitions has to be now written as 



T^cS — ^ C'xOx 



X 

where Cx are the Wilson coefficients computed in the MSSM and Ox are the same four-quark operators 



as in eq. (2.1). Assuming for definiteness that sparticle masses are of the order of AfsuSY = 500 GeV, 
and taking into account the QCD RG running of Cx between MgusY and the hadronic scale as well 
as matrix elements of the operators Ox between the meson states in the manner described in [|2^ one 
obtains for the supersymmetric contribution to the K^-K'^ transition amplitude: 

{K^\n.^\K') « M^Fl [0.15 (Cs^u^s^ + C^^^y) " 6.0 (Cf^J^Y + C|^s\) 

-11.5 (Cs^^-sY + Cs™sy) - 13.84 C^^f^ + 22.48 C|{;«y] (5.7) 

where we have used as{Mz) = 0.1185. The large numerical factors]^ in the second line originate from the 
RG running and from the chiral enhancement factor (Mif/(ms -t-m^))^ « 18 for ms(2 GeV) = 110 MeV. 



For the supersymmetric contribution to the B^-B^ transition amplitude one has to replace in eq. (5.7) 
MkF^- by Mb^F^^ (g = d or s) and the numbers in the square bracket by: 0.24, -0.49, -0.94, -0.97 and 
1.27, respectively. Note, that there is no chiral enhancement in this case as (Ms,/(mb -I- m^))^ « 1.65. 
Using the standard formulae 



AMk = 2 Re {K°\nes\K") 

i7r/4 

V2AMk 
AAf, = 2|(B°|Heff|i?°)l 



and plugging in numbers one finds 



AMk = 3.87 x 10'' Re x [. . .]) ( ) ps 



-1 



Ml 



® Except for sneutrinos whose mass squared matrix consists of the LL 3x3 block only. 

®Some uncertainties of order few percent in these numbers are due to the uncertainties of the factors parameterizing 
matrix elements of the operators Ox for X=SLL, SRR, VLR, SLR, TL, TR. ^ 
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Figure 8: Contribution of gluino-squark box diagrams neutral meson mixing. Crosses denote mass 
insertions. 

SK = -2.58 X 10*^ Im (m| x [. . .]) f^^^) e'^'^ (5.8) 



AMd = 6.45 X 10^ Re (Af| x [. . .]) ( ) ps 



-1 



Ml 



2 



AM, = 8.78 X 10^ Re (M| x [. . .]) ( ) ps" 



where M| is some average mass of squarks and [. . .] denote the content of square brackets from eq. (5.7) 



appropriate for a given transition. 

In the lowest order in the mass insertion approximation each of the Wilson coefficients Cx for | AF| = 2 
transitions like K°-K° and B°-B° can be represented as a product of a function of M| and rrig and of 
two mass insertions defined as ^ |l5j 

\<^xy) 

where X,Y — L,R, and the indices J, I label generations. Neglecting contributions other than those 
generated by gluino exchanges, one has 

Clhh « 0.1 [a' {5^,,y"{5U" + b' {5^lh)"{5^hl)"] (5-9) 

etc. [|0l, where J I = 21, 31 and 32 for K°-K°, B°-B'^ and B°-B° transitions, respectively. 

Comparison of the numbers in eqs. ( |5.8| ) with the experimental values: AMk = 0.0053/ps, ex — 
2.28 X 10"'^ and AM^ = 0.496/ps illustrates the so-called supersymmetric flavour and CP problem: taking 
into account that the dimensionless factors M~ x [. . .] in eqs. ( |5.8| ) are 0{1), it is clear that the typical 
contribution to AMk, £k, AMb^ and to many other measured quantities like e'/e, BR{B Xs'j) 
pl| , [Tsl , etc. in the MSSM with the flavour and CP violation in squark mass matrices is several orders 
of magnitude too big. Any theory of supersymmetry breaking has to face the problem of explaining the 
smallness of the mass insertions S'^y- 

Adopting the rough criterion that the gluino contribution alone to |AMif | and {skI should not exceed 
the experimental values of these quantities (and barring possible cancellation between different mass 
insertions) one obtains for small and moderate values of tan (3 the limits shown in the middle column of 
table 0. For comparison in the first column we show the limits obtained in the paper [ pO[ . The differences 
stem from slightly different treatment of the NLO QCD RG evolution and of the matrix elements of the 
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operators involved (our approach is based on ref. p2[ ) but are inessential for the order of magnitude 
estimates of the limits. 



Table 1: Upper limits on mass insertions obtained from ex for Mq = 500 GeV. The limits scale approx- 
imately as M|. X = {rrig/Mq)^. Limits on [((5f^)^^] are the same as for [(^i2)iL]- As follows from 
numbers in the first two of eqs. (5.8), the corresponding limits on real parts of the product of insertions 
are simply 12.5 times weaker than those given below. (This simple rule is not satisfied by the numbers 
quoted in ref. M.) 





ref. |0) 
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g 






low tan /5 


tan /3 = 50 


X 




1 


0.3 


2.9 X 10^3 


2.7 X 10-3 


2.5 X 10-3 


1.0 


6.1 X 10^3 


6.0 X 10-3 


1.7 X 10-3 


4.0 


1.4 X 10^2 


1.5 X 10-2 


2.0 X 10-3 


9.0 




1.4 X 10-2 


2.4 X 10-3 


X 


Vlim [('5f2)i J 


1 {m)LR\:»\{5?2)RL\) 


0.3 


3.4 X 10-* 


2.5 X 10-4 


2.2 X 10-4 


1.0 


3.7 X 10-4 


2.9 X 10-4 


2.2 X 10-4 


4.0 


5.2 X 10-4 


4.2 X 10-4 


2.5 X 10-4 


4.0 




6.5 X 10-4 


6.9 X 10-4 


X 


^|Im 




1 


0.3 


1.1 X 10-4 


8. 


2 X 10-^ 


8.0 X 10-^ 


1.0 


1.3 X 10-4 


9. 


5 


X 10-^ 


9.2 X 10"^ 


4.0 


1.8 X 10-4 


1.4 X 10-4 


1.3 X 10-4 


9.0 




1.9 X 10-4 


1.8 X 10-4 



It is interesting to note, that because Sk puts stringent bounds only on imaginary parts of products 
of two mass insertions, bounds on almost real and almost imaginary mass insertions are provided only by 
AMk and are order of magnitude weaker, although such a conspiracy seems not very natural. Stronger 
absolute bound on the imaginary part of the mass insertion itself exist only for (S^j^)^'^ and follows from 
s'/s: |Im(<5f^)i2| < 10-5 QQ. 

For mass insertions generating transitions between the third and the first two generations of quarks 
only much weaker bounds are available. Limits on (Sxy)^^ insertions from the giuino box contribution to 
AMd have been derived recently in ref. |Q. Similar limits on (i^xy)^^ insertions will become available 
once AMs is measured. At present however, stringent bounds from the B — > decay exist only for the 
insertions (5^^)23 and {Sj^j^)'^^: < 0.07 x (Mq/l TeV). The remaining insertions are bounded 

rather weakly |2|. 

For large tan/3 the standard analysis of bounds on mass insertions derived from |Ai^| — 2 transitions 
based on giuino box diagrams of figure || is not sufhcient. Scalar flavour changing neutral Higgs boson 
couplings can generate additional contributions oc tan4 fj to the Wilson coefficients C^^^ through the 
double penguin diagrams of figure || (contributions to C^^^ , and C^^^ Wilson coefficients are suppressed 
because of the mutual cancellation of and A" contributions) and these contributions have to be taken 
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into account. Dominant source of the flavour changing neutral Higgs boson couplings in the case of 
flavour violation in squark mass matrices are the diagrams b) and c) shown in figure ^. Calculating those 
diagrams one gets the couplings 



(5.10) 



where the matrix coefficients X^j^ — 



i^LRy are given by 



s 

RL 



JI 



x'^ tan^ (3 



eas 



xD{mlMh 



1.1 rridj 



Mw Mw 



u 



^RR) 



cos a, — sin a and isin/3 for = iJ", or A", respectively and D{a,b) some dimensionless 



with 
function. 

It turns out that even the limits on {SEl)^^ and (^/j^)^^ mass insertions are affected by the double 
penguin contribution which is significantly enhanced by the big numerical factor multiplying CgyJ^ in 
eq. (5/7). The effect of double penguin contribution is seen in table |l| in the limits on imaginary (and 
real) parts of [((Jf^,)^^]^ and [((^/j/j)^^]^ which for rrig > Mq become stronger by one order of magnitude 
compared to similar limits for lower tan (3 values. That the improvement is seen only for rrig > Mq follows 
from the fact that the couplings ( 5.10| ) are proportional to rUg. The limits on (^Li)^^(5ijfl)^^ are not 



improved because, as is clear from eq. (5.9), the gluino box contribution to Cgygy contains already a 
term proportional to (<5iL)^^(^i{fl)^^- 



In the same manner, bounds on the insertions {S^l)^^ and {SEr}^^ (and when AMg is measured 



also on (5fi)^^ and {S^j^)^^) derived from B°-B° mixing Q should also be modified for large values of 
tan/3. We have found, however, that the actual bounds depend also the chargino box and double penguin 
contributions and can not be therefore presented in a simple way. 

Another interesting effect related to the flavour changing couplings ( 5.10| ) generated for large tan/3 
by non-zero LL and/or RR mass insertions is a growing like tan^ /3 contribution to the amplitudes of 

~^ f^^l^~ and/or fJ-^ fJ-~ decays [|l^. Calculating the contribution of the diagram shown in 

figure H with the couplings ( 5.1C| ) one finds for the coefficients a in the amplitude (4.5) 



^ tan^ (3 



xZ?(m? M?) 



mg^. 
Ml 



RR) 



where q = d = 1 and q 



s = 2 for _BJ] and decays, respectively, and the coefficient b in the amplitude 
(4.5) is given by the similar expression with + changed to — in the square bracket. It has been shown 
p9[ , that for tan/3 ~ 50, Ma ~ 200 GeV and with mass insertions of order 0.1 the branching ratios 
predicted in the MSSM can exceed by one or two orders of magnitude the present experimental limits 



(4.6). It is however important to check whether this remains true when all the available constraints are 
respected, including the ones imposed by AMs and AMd and taking consistently into account the double 
penguin contributions to these quantities. The results of this exercise are shown in figure ^ where we 
show the branching ratios of the decays 5° l^'^IJ'~ and B^ M"*"/^" versus the mass difference AMg 
for a sample of the MSSM parameters varying the insertion (SEl)^^ in the range (0.01,0.1). All points 
giving rise to experimentally unacceptable values of AM^ and/or BR{B Xgj) have been discarded 



Points for which < have been discarded as well. Since we set the insertions and ((5 



^RR. 



32 
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Figure 9: Correlation of BR{B^ (left panel) and BR{B'^ IJ'^I^') (right panel) with AMg 

in the MSSM with flavour violation in the squark sector. The single nonzero mass insertion ((5^^)^^ has 
been varied in the range (0.01, 0.1). tan/J = 50 and Ma = 200 GeV. 
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equal to zero the increase of BR{B^ — > ji^ compared to the SM prediction seen in the left panel of 
figure ^ is mainly due to the effects discussed in section 4. 

It is clear form figure I that even with all the cuts imposed the possible values of the branching 
ratio BR{B^ — > can still be above the present experimental limit ( |4.6| ) which means that the 

non-observation of the B^ — > /x^/i^ decay imposes nontrivial constraints on the MSSM parameter space 
and on the mass insertions (Sj^j^)^^ and (S^j^)^^ ■ 

Finally, comparison the possible effects in the MSSM without and with flavour violation in the squark 
mass matrices (figures ^ and respectively) leads to the interesting conclusion that, within the super- 
symmetric framework, observation of the B'^ decay at the level close to the present BaBar limit 



(4.6) (i.e. with BR above ~ 6 x 10~^ in general and above ~ 3 x 10~^ if unnaturally large and very 
unlikely values of the stop mixing parameters At are not taken into account), apart for implying that 
the scale of the Higgs sector is not far from the electroweak scale, would be a very strong evidence of 
non-minimal flavour violation in the quark sector. 

6 Effects of flavour violation in the lepton sector 

To complete the picture of flavour violation in the supersymmetric extension of the SM model we discuss 
briefly also the lepton sector. 

To account for the observed atmospheric and solar neutrino oscillations |Q the analog of the CKM 
mixing matrix, the so-called Maki-Nakagawa-Sakata mixing matrix Jssf , has to be introduced in the 
leptonic sector of the SM or the MSSM. Under the assumption that the mixing occurs only between the 
three know neutrino flavours (no sterile neutrinos) which is supported by the SNO results , the MNS 
matrix U is of dimension 3x3 and is usually parameterized in a similar way as the CKM matrix 

(C12C13 S12C13 Si3 \ 
S23C13 (6.11) 
C23C13 / 

where C12 = cos ^12 etc. and where we show only the entries directly related to observed oscillations 
and neglected all possible CP violating phases. Non-zero angles 9i2 and ^23 are responsible for solar and 
atmospheric neutrino oscillations, respectively. 

The pattern of U emerging from the experimental data ^ 

\Uii\ « IC/12I «l/%/2 

IC/23I « IC/33I ~ I/V2 (6.12) 
IC/13I «0 

is called bi-maximal mixing and is distinctly different from the hierarchical pattern of the CKM matrix. 



Nontrivial mixing matrix U (6.11) induces of course also flavour violating processes with charged 
leptons, such a /i — > 67 or — > e/i etc. but at rates which are unmeasurably small (e.g. BR(^ —> 67) < 
-[^g-50-j jtqj, neutrino mass squared differences required to explain the Superkamiokande data: 

Aml^ = ml^ - ml^ « 3.2 x IQ-^ eV^ , 
Amli = ml^ -ml^^O (IQ-^) eV^ 

and masses compatible with constraints imposed by cosmology (X^a'^i'a < ^w eV). Thus, if the MNS 
mixing matrix is the only source of flavour violation in the leptonic sector, neutrino oscillations remain 
the only observable lepton flavour violating phenomenon. 
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In supersymmetric extension of the SM lepton flavour violation can originate also in the slepton mass 
squared matrices. Existing experimental upper bounds: BR{^ —> 67) < 10""'^^, BRir — > e(/i)7) < 
10"^ put stringent constraint only on the mass insertion | ((5|^fl)^^ | which has to be smaller than 10^^; 
constraints on the other sleptonic mass insertions are of order fewxlO^^ p^ . 

Interesting links between the lepton flavour violation originating in the slepton sector and neutrino 
masses and mixing exist in the see-saw scenario in which observed small neutrino masses result from 
exchanges of right-handed neutrinos vr of masses Af^^ ~ 10^" — 10^'' GeV in the GUT-type framework. 
Firstly, the RG running of the parameters of the theory between the scales Mqut and M^^ necessarily 
induces lepton flavour violating mass insertions. It turns out that the experimental limits on /i — s- 67, 
T — s- e(/i)7 decays put interesting constraints on realizations of the see-saw mechanism in the GUT-type 
scenarios pO| . Secondly, lepton flavour violating originating in the slepton sector can influence neutrino 
masses and mixing via quantum corrections below the scale M^^. Let us discuss this point in some 
details. 

Quantum corrections to neutrino masses and mixing below the M^^ scale are of two types. The 
first one are the corrections depending on \d.{Mi,j^/Mv/) which are accounted for by integrating the 
renormalization group equations pl| of the dimension 5 operator between the scales Mi,^ and Mw- The 
most interesting aspect of the RG equations is their fixed point structure : whenever the RG running 
is substantial, the mixing angles evolve in such a way that at the scale either [/31 = or C/32 = 0. 
In both cases one gets the following relation between the mixing angles 

^i^' 2012 = .2 J\ . ^2 Sin' 2023 
Because of the CHOOZ limit 5^3 < 0.16 |39] this is incompatible with the bi-maximal mixing pattern 



( 3.12 ) favoured by the solar and atmospheric neutrino data. This means that always for exact three-fold 
or two-fold degeneracy of neutrino masses at the scale M^^^ , or for approximate degeneracies of neutrinos 
having the same CP parities, when the RG running is substantial |Q, the mixing angles obtained from 
the see-saw mechanism are phenomenologically unacceptable. Note also that these are precisely the most 
interesting cases: three-fold degeneracy of neutrinos will be required if the neutrinoless double beta decay 
is found at the level requiring m^*^ 0.5 eV. More generally, see-saw scenarios giving naturally large 
mixing angles may be easier to find if the spectrum of neutrino masses is (approximately) degenerate. 

The unacceptable pattern of mixing generated by RG running can however be changed by the second 
type of quantum corrections to the neutrino mass matrix - the so-called low energy threshold corrections 
- if there is some lepton flavour violation in the slepton sector p3[ . 

In the basis in which the neutrino mass matrix (m^)'^^ generated by the underlying see-saw mecha- 
nism is diagonal the corrected neutrino mass matrix can be written as ^ 

m^r^ + [U^ (l^m^^^ + m(0) j) u] (6.13) 
where U is the uncorrected MNS matrix and 

jAB ^jAB _gABjA^ (6.14) 

summarizes the RG [I,^] and low energy threshold (/^j^ ) corrections. The most interesting part of the 
latter corrections take the form 

Itr-('5U^''x/(Mf,m?,) 
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where the function of chargino and slepton masses / is typically of order fewx (10 ^ — 10 (contributions 

r 

rg 



of (5^^ and S^i^fj to are smaller). For comparison, for My^ ~ 10^° GeV, « 10 '"^ x ta.T? [3, and 



/|"g, /;;g are negligible. 

In the case of the (approximate) degeneracy of the zeroth order neutrino masses rn'y'} w in'^^ the 
matrix U is fixed by the condition 

Y,UAalABUBb = Q (6.15) 

AB 

(the freedom U ^ U ■ Raby where Rab is an arbitrary rotation of the Va and z^f, neutrino fields, is used to 
diagonalize the "perturbation"). This leads to the fixed point-like relations between the mixing angles 
which are different than the RG evolution provided ^ \lfgc\ @]- 

As an example consider initial degeneracy of the three neutrinos m^^ « m^i^ ~ —"ni^^ and only one 
dominant correction I'^ . In this case interesting results are obtained for m,yj « m^., (or rrii,^ « n^'j^s) 



and dominant I^^ correction (i.e. (S^^i^)'^^ ^ 0). The condition ( 6.15 ) then gives 



Si3 = - cot 26*23 tan (cot 6*12) 

which is compatible with the bi-maximal mixing and small C/13 — S13 element. Moreover, for the mass 
squared differences one obtains 



Amgoi = -Ami cos 26*12 sin 26*23/: 



th 



th 



Aml^^ = -Amlil + cos^ 6*12) sin 26*23/: 
that is, for the bi-maximal mixing: 

Am32 > Am2i ^ 

in agreement with the experimental information. Am^^^-^ sa 3 x 10^^ eV^ requires then mj/ « 1 eV^ and 
(^ll)^"^ 0.5 with the interesting implications for the t ^ iij decay. Amg^j of right magnitude can be 
generated either by departure from ^12 — t:/A or, by another, hierarchically smaller, correction: flavour 
conserving I^^ — 1^6^^ with either /^ 7^ or /^ 7^ (e.g. /^ 7^ from RG running for not too large 
value of tan/3), or flavour violating correction /j^, or /^^. 

There can be, of course, other interesting cases with more complicated interplay of RG and low energy 
threshold corrections 



7 Summary 

We have reviewed recent developments in exploring flavour dynamics in the supersymmetric extension of 
the Standard Model. Emphasis has been put on possible interesting effects in 6-physics arising for large 
values of tan (3 and not too high a scale of the MSSM Higgs boson sector, both in the case of minimal 
flavour violation and in the case of flavour violation originating in the sfcrmion sector. We have discussed 
the importance of the flavour changing neutral Higgs boson couplings generated by the scalar penguin 
diagrams and their role in constraining the MSSM parameter space. We have shown that in the case 
of minimal flavour violation the experimental lower limit on B^-B^ mass difference constrains branching 
fractions of the decays 5° ^ — > ^^/^^ possible in the MSSM. We have also pointed out that observation 
of the B^ /^^M^ decay with BR at the level 3 x 10~* (and even lower if AM^ turns out to be 
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bigger than 15/ps) would be a strong indication of nonniinimal flavour violation in the quark sector. 
Flavour violation connected with neutrino oscillations has been also discussed. It has been argued that 
in some physically interesting situations flavour violation originating in the slepton mass matrices can be 
responsible (at least in part) for observed pattern of the neutrino mixing and mass squared differences. 
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